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IMPLEMENTATION OF INTERACTION ALGORITHM TO NON-MATCHING 
DISCRETE INTERFACES BETWEEN STRUCTURE AND FLUID MESH 


PO-SHU CHEN* 


Abstract. This paper presents software for solving the non-conforming fluid structure interfaces in 
aeroelastic simulation. It reviews the algorithm of interpolation and integration, highlights the flexibility 
and the user-friendly feature that allows the user to select the existing structure and fluid package, like 
NASTRAN and CLF3D, to perform the simulation. The presented software is validated by computing the 
High Speed Civil Transport model. 
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1. Introduction. 

1.1. Background. The importance of aeroelastic problems has been widely recognized in many en- 
gineering fields like acoustics problems, airfoil oscillations, and flutter predictions. Since the aeroelastic 
analysis considers not only the properties of fluid but also the flexibility of the structures, it improves the ca- 
pability for designers / analysts to understand the interaction of fluid/structure, which improves the accuracy 
of preliminary and design loads and leads to a reduction in development and production costs. 

However, the analysis of aeroelasticity involves solving fluid and structural equations simultaneously. Be- 
cause most aerospace vehicles are often dominated by large structural deformations, fully coupled procedures 
are required for accurate simulations. 

Different methodologies have been developed for computational analysis. The first class is t ightly coupled 
aeroelastic analysis, i.e., solving both structures and fluids in a single computational domain. The major 
disadvantage of this methodology is the ill-conditioned matrices associated with two physical domains. The 
secondary disadvantage is not being able to use the existing CFD codes. There has been a large investment 
of time and money in the development of classical, rigid CFD programs that have been t ailored specifically to 
different applications. A tightly coupled procedure is not able to take the hill advantage for these specialized 
and well-trusted programs. 

On the other hand, the loosely-coupled methodology uses two independent disciplines by exchanging data 
at interfaces between fluids and structures. This allows it to take full advantage of existing, well-developed 
programs like NASTRAN for structure analysis and CFL3D for fluid analysis. A completely aeroelastic 
simulation cycle could be described as in figure 1.1 and a typical simulation may need about three to five 
cycles. 

Obviously, two different disciplines will have non-matching discrete meshes due to their different interests. 
For example, the fluid mesh may have a finer grid at the wing tip to catch the phenomenon of vortex, while 
the structure grid lias a relatively coarse grid since the wing tip is not the area of stress/strain concern. 
Several approaches have been proposed in the past for solving the fluid /structure interaction problems on 
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Fig. 1.1. Typical Aeroelastic Simulation Cycle 


moving and deforming meshes. 

The motivation to develop a package is to improve the aeroelastic simulation conduct by the Multi- 
disciplinary Optimization Branch at NASA Langley Research Center. FASIT, which stands for fluids and 
structures interface toolkit, developed by Prof. Marilyn Smith, is currently used for interpolation and in- 
tegration between fluid and structure analysis. However, this code is generally difficult to use. The other 
disadvantage is the geometry definition, which prevents the accurate calculation for any object but the wing. 
The new program, LMT. has been developed to be a “bridge” between CFD and FEM software for aeroe- 
lastic simulation. LMT stands for Load and Motion Transfer program. It is able to interpolate the initial 
nodal coordinates of the fluid mesh from the structure nodal displacement, and to integrate the structure 
nodal force from the fluid pressure. The algorithm behind this program was proposed by Prof. Charbel 
Farhat and Michel Lesoinne at University of Colorado, Boulder. 

1.2. Goals. The design and implementation of this new package are guided by several principles. These 
goals are described as follows. 

1.2.1. User Friendliness. The new package has to be easy to use and straightforward, with no need 
to convert data to different formats and no need to specify geometry reference 1 points. 
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1.2.2. Flexibility. Allows user to take or to switch different CFD/FEM packages easily, thus re- 
searchers are able to select the most appropriate software for loosely-coupled aeroelastic simulations. To 
achieve this goal, the code must, be able to understand, at least,, major CFD formats like Plot3D or Tot Plot . 

1.2.3. Extensibility. Extensibility allows the program to he equipped with the latest integration 
method, or different data format for a new CFD/FEM program, with only minor modification of the code. 

1.2.4. Accuracy. The algorithm enforces the satisfaction of conservation of momentum and energy. 

2. Algorithm. To ensure the quality of the transfer, a good algorithm has to preserve the consistency 
and conservation. The consistency requires that the summation of the nodal force vector on the struc- 
ture mesh must be equal to the resultant force and moments induced by fluid pressure on the fluid mesh. 
Conservation refers to the virtual work performed by the load vector on the structural mesh with virtual 
displacement equal to the work performed by fluid pressure on the fluid mesh with the associated virtual 
displacement. 

A brief review of the algorithm is presented here. The first section is the load transfer algorithm while 
the second section is the motion transfer algorithm. 

2.1. Load Transfer Algorithm. Let u refer to the admissible virtual displacement function. Subscript 
F refers to the fluid domain while S refers to the structure domain. I denotes the interface between structure 
and fluid domains. The trace of up and its satisfy 


(2.1) u f = its on P. 

We could describe the displacement of every surface point in the fluid mesh as a function of the nodal 
displacements of the structure model as follows. 

i=is 

( 2 . 2 ) up, = YL C ‘> "*■ J e Tf ' * € r - s ’’ 

*= 1 

U F is the discrete value ofu F . Similarly, its, is the discrete value of Us- C u are constants which depend on 
the approximation method. 

The virtual fluid displacement function is discretized as follows: 


J-JF 


(2.3) 


it l = D J “P J € r 

j= 1 


F- 


The virtual work on l> by the action of the fluid pressure force is 


(2.4) 

(2.5) 

(2.G) 


SW F = / 

J r, 


-]m) up ds 


J=JF r. 

= V / (-pn)DjU F jds 
j = ' ir - 
j=j( 
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3>j has the physical meaning of numerical pressure flux. 


(2.7) 





D , ds 


The virtual work on Ts- by the action of the structure force could be written as 


( 2 . 8 ) 


>-is 

^W’s = E “Si • 

i=l 


To satisfy the principle of energy conservation, JIT> = fll's-, we conclude that 


(2 - 9) fi=E^ c J‘- 

j= i 

The first term, pressuie flux, is independent of the structure code, while the second term depends only 
on the approximation method. 

Since the finite element method has dominated the solution method of the structure problems, the 
structural element displacement field on Ts is expressed as 


( 2 . 10 ) 



USi - 


Combine Eq.(2.10) with Eq.(2.1), we have 


(2- 11 ) U Fj - UF(Sj) - Us(Xi) = E N i(Xj) j G I f , i G Ts. 

i= 1 

Following Nj(\j) = Cjj , Eq.(2.11) could bo expressed as 

( 2 - 12) Si = E 

7=1 

This is the formula adopted in the new package. To compute the Ni(\j), not only the structure nodal 
coordinates but also the structure element topology have to be provided. 

2.2. Motion Transfer Algorithm. To transfer the motion from structure to fluid surface, recall Eq. 

( 2 . 2 ), 

i=is 

u Fj - ^2 US > j € T/r, i £ T.s. 

t= 1 

Similarly, if we choose the shape functions for the approximation as the load transfer, the above equations 
become 

•=*s 

( 214 ) «f> = E N ‘j Us > ■> e O - i G T s . 

/=i 
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Fluid Mesh 


Structure Mesh 


3. Implementation. The procedures of load/motion transfer are divided into two parts: projection 
pha.se and integration/ interpolation phase, with one program associated with each phase. 

3.1. Projection. The projection code, matcher, project fluid nodes to the structure surface, and then 
computes the position of the projected fluid nodes on the associate structure element in terms of natural 
coordinates. This program needs to be executed only once at the beginning of simulation as long as the 
deformation is fairly small or for the problems with different, surface load or deformation but the same mesh 
definition. 

3.2. Integration and Interpolation. The second program, LMT, takes the structure/fluid coordi- 
n;in-. si incline displacement or fluid pressure, and the natural coordinates file created by matcher, to do 
tin intrgiation or interpolation. 

3.3. Norm problem. One of the tricky parts regarding implementation is the direction of the norm, 
i ni example, on the upper surface of the wing, we desire the downward norm since the pressure force is 
dow nu.ivil too. On the other hand, we desire the direction of the pressure force upward on the lower surface 
sun . i i m lnwer surface provides the lift. 

In case of an unstructured fluid mesh, the problem is trivial. Since the boundary condition has to 
hr explicitly given, we can arrange the boundary facade counter clock wise as seen from the inside of the 
structure, then the norm vector can he computed accordingly. 

For structured mesh, however, the boundary facade is implicitly given. The user may not. even know the 
node number but the indices of the mesh. Therefore, there is no difference of the upper surface and lower 
surface from the numerical point of view. A special flip option is implemented to indicate whether the norm 
vector for each zone needs to be “flipped” or not. 

4. Some Numerical Results. The capability of the program has been demonstrated by solving the 
high-speed civil transport (HSCT) model. The answers are verified by FASIT. 

4.1. HSCT model. For the structure model, the number of nodes is 226 and the number of triangular 
shell elements is 1274. The fluid mesh is structured with four zones surrounding the structure model. These 
four zones are upper/lower wing and upper/lower fuselage. Figure 4.1 shows the fluid mesh. 




Fig. 4.1. Structure and Fluid Grid 


4.1.1. Load Transfer. The results show good consistency with the FASIT code at the upper and lower 
wing with only one percent of difference in the z direction. The force on upper /lower fuselage can not be 
verified by FASIT due to the geometry definition limitation. However, these forces cancel each other out as 
we expected. 


LMT Result FASIT Result 


Zone 

upper wing 
upper fuselage 
lower wing 
upper fuselage 


-3.8639e+05 lb 
-1.4642e+05 lb 
5.9174e+05 lb 
1.4236e+05 lb 


-3.8653e+05 lb 
N/A 

5.9119e+05 lb 
N/A 


4.1.2. Motion Transfer. As in figure 4.2, the smooth deflection along the wing root and fuselage 


5. Discussion and Conclusion. LMT provides an ideal tool for aeroelastic simulation. It could serve 
as the testbed for different integration methods, or as the tool for people who need to have a quick answer 
for aeroelastic problems. 
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dz 

2.500000 
2.307143 
2.1 14286 
1.921429 
1.728571 
1.535714 
1.342857 
1.150000 
0.957143 
0.764286 
0.571429 
0378571 
0185714 
- 0.007143 
- 0.200000 



Fig. 4.2. Deflection Comparison 


Some future improvements include 

• Viscous force. 

The structure nodal force induced byfluid is composed of two parts, i.e., pressure and viscous force. 
This package considers the pressure force only. Viscous force will be added in a future release. 

• Unstructured Fluid Mesh. 

Most of the existing CFD packages use structured fluid meshes. However, unstructured fluid mesh's 
are gaining popularity these days due to their less stringent memory requirement, and greater flexi- 
bility for the area of interest. The future release of the LMT package will allow unstructured fluid 
meshes. 

• Two-dimensional Problems. 

The program is designed for three-dimensional aeroelastic simulation. However, it will be expanded 
for two-dimensional aerolelastic problems also. 

• Different Integration Methods. 

The reasons to choose this algorithm are accuracy and simplicity. However, other algorithms, like 
finite-plate spline, biharmonic-inult.iquadric method, could be added easily for research purpose's. 




• Different Data Formats. 

LMT recognizes PLOT3D and NASTRAN formats only. Other different data formats are desired 
to increase the flexibility of the code. 

• Different Elements. 

The only stiuctuie boundary facade allowed is triangular at this moment. However, complex prob- 
lems involve a large variation of different elements. Quadrilateral, beam, and other type of elements 
will be added soon. 
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